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I. INTRODUCTION 

In the present paper, a particular approach to understanding lepton mixing, proposed in and furthcr studied in 
0, will be examined in more detail. First, we briefly review the approach. 

Of course, the Standard model interaction term for /3 decay or n" —f erVe includes the leptonic piece: 



C = -^W^p eLli^Ve + h.c, 
The object is now known to be a linear combination of neutrino mass eigenstates, /5^: 

where, in a basis with the charged leptons diagonal, the fuU lepton mixing matrix is written as: 

Kel Ke2 Ke3 

K = I Kf,i K^2 Kf,3 

Krl Kr2 Kr3 



(1) 
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(3) 



As has been discussed by many authors [iT|-[22| the results of neutrino oscillation experiments are (neglecting possible 
phases to be discussed later) consistent with the "tribimaximal mixing" matrix: 
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Many different approaches have been used to explain the form of K. A "natural" , and often investigated one uses the 
parallel three generation structure of the fundamental fermion families as a starting point. An underlying discrete 
symmetry S3, the permutation group on three objects, is then assumed. [l^l-jl^] The permutation matrices S are, 
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This defining representation is not irreducible. The 3-dimensional space breaks up into irreducible 2-dimensonal and 
1-dimensional spaces. One may note that the tribimaximal matrix, K^bm is an example of the transformation which 
relates the given basis to the irreducible one. This faet provides our motivation for further investigating the 
symmetry, even though many other interesting approaches exist. Of course, the symmetry requirement reads, 



[S, M,] = O, 

where S stands for any of the six matrices in Eq.((5|) and My is the neutrino mass matrix. 
By expUcitly evaluating the commutators one obtains the solution: 



= al + iid. 
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a and /? are, in general, complex numbers for the case of Majorana neutrinos while d is usually called the "democratic" 
matrix. 

It is easy to verify that this My may be brought to diagonal (but not necessarily real) form by the real orthogonal 
niatrix, R = Ktbm defined above: 



R^{al+f3d)R = 

R may be written in ternis of the eigenvectors of M^, as: 

R=[ri 



a O O 
O a + 3/3 O 
O O a 



r2 rs 



(8) 
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For example, fi is the first column of the tribimaximal matrix, Eq.(|3]). Physically one can assign different masses 
to the mass eigenstate f2 in the 1-dimensional basis and to the (doubly degenerate) eigenstates fi and in the 
2-dimensional basis. At first glance this sounds ideal since it is well known that the three neutrino masses are grouped 
into two almost degenerate ones ("solar neutrinos") and one singlet, with different values. However, since we are 
demanding that R be taken as the tribimaximal form, the physical Identification requires fi and r2 to be the " solar" 
neutrino eigenstates rather than the degenerate ones ri and f^. This had been considered a serious objection to the 
present approach since often a scenario is pictured in which the mass eigenvalue for riis considerably larger than the 
roughly degenerate masses associated with fl and r^. A way out was suggested in [l[ where it was noted that, for 
values of mi +7112 + ^13 larger than around 0.3 eV, the neutrino spectrum would actually be approximately degenerate. 
This may be seen in detail by Consulting the chart in Table 1 of [l| wherein the neutrino masses are tabulated as a 
function of an assumed value of the third neutrino mass, m^. Actually it is seen that there is also a region around 
TO3 ~ 0.04 eV and toi +7712 + 7713 ~ O.lSeV where an assumed initial degeneracy may be reasonable. To make physical 
sense out of such a scenario, it was suggested that the neutrino mass niatrix be written as. 



(10) 



where Mjf'' has the fuU S3 invariance and has degenerate (at least approximately) eigenvalues. Furthermore, the 



smaller M^^-* is invariant under a particular 5*2 subgroup of 5*3 and breaks the degeneracy. Finally, M^^' is invariant 
under a different ^2 subgroup of 5*3 and is assumed to be smaller still. The strengths are summarized as: 

Mi°'> > M^i) > Ml'^l (U) 

This is inspired by the pre-QCD flavor perturbation theory of the strong interaction which works quite well. In that 
case the initially unknown strong interaction Hamiltonian is expanded as 



(2) . 



(12) 



Here iJ*^"^ is the dominant SU{3) flavor invariant piece, H^^'> is the smaller Gell-Mann Okubo perturbation [2^ 
which transforms as the eighth component of a flavor octet representation and breaks the symmetry to SU(2) and 
H^'^\ which transforms as a different component of the octet representation and breaks the symmetry further to the 
hypercharge U(l), is smaller still. 

There is a possible immediate objection to the assumption that the neutrino mass eigenvalues be degenerate in the 
initial S3 invariant approximation; after all Eq.(I5]) shows that there are two different eigenvalues a and a + 3/3. This 
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was overcome by recognizing that these are both complex numbers and that they could both have the same magnitude 
but difFerent directions. Having the same magnitude guarantees that all three physical masses will be the same. This 
introduces a physical phase corresponding to the angle between a and a + 3/3. 

In the strong interaction case, the initial SU(3) invariance was found to be reasonably weU obeyed. It is thus natural 
to ask what prcdictions may exist in the initial S3 invariant approximation in our neutrino model. It was found [l| 
that the leptonic factor for neutrinoless double beta decay, rriee could be predicted in this limit to be, 

|mee| = y \/5 + 4:C0S1p, (13) 

where m is the degenerate neutrino mass and i/j is the Majorana type phase mentioned above. This led to the 
inequality 

m > |mee| > m/3. (14) 

The next step in the program is to consider the efFect of the perturbation M^^-*. Many authors [2^-[3l[ have 
suggested that a /i-r symmetry ((23) symmetry in the present language) is associated with tribimaximal mixing in 

the neutrino sector. Thus it is a natural iS'2 symmetry choice for Recently, Chen and Wolfenstein applied 

this type of perturbation to our present model with the additional assumption that the Majorana phase ip takes the 
value TT. This corresponds to CP conservation. Their result for \mee\ is in agreement with the lower limit in Ea. (jl4p . 
Here we will investigate the first perturbed case without assuming that special value of ■0. 

Before going on to this we will present an amusing argument to show that the (23) perturbation is naturally 
associated with the tribimaximal form (modulo the majorana type phase ^/;) rather than a tribimaximal form multiplicd 
by a rotation in the two dimensional degenerate subspace (which is physically irrelevant at the invariant level). 
This is based on the fact that degenerate perturbation theory must be employed, which leads to a stability condition. 
Further we will show that other S2 perturbations are mathematically consistent but do not lead to the desired 
tribimaximal form. 



II. EFFECTS OF DIFFERENT PERTURBATIONS 



In the present framework there are three different possible perturbations, each characterized by the S2 subgroup 
which remains invariant. Let us first consider the favored perturbation which leaves invariant the S2 subgroup, 
consisting of S*^^-* and S'^^'^^. Apart from a piece which may be reabsorbed in Eq.([7|), such a perturbation has the 
form, 



(15) 



where t and u are parameters. It is convenient to adopt the language of ordinary quantum mechanics perturbation 
theory. We should then work in a basis like Eq.([8]) where M^, in Eq.([7|) is diagonal. However, because of the double 
degeneracy between the eigenvectors ri and fa in Eq.(|3]), the matrix R is not the unique one which diagonalizes M,^. 
We should really use the more general matrix RX[S^) where X{S^) is given by: 
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In this basis A has the form: 
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(17) 



+ c^{t ~ u) 



Here, c = cos^ and s = sin^. Note that, before adding a perturbation, the 6*3 symmetry predicts the lepton mixing 
matrix to be RX{(^) rather than the desired tribimaximal form, R. 

In perturbation theory, the first correction to the m*'* eigenvector involves the ratio ^"J,'^_ J,'"^ . For degenerate 



perturbation theory it is of course necessary that the numerator vanishes for those states with E', 
simply require for the (13) matrix element: 
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This yields in general, sin{2£) = 0. The solution with ^ = O is the desired tribimaximal form. The solution with ^ = tt 
just changes the signs of the first and third columns. However, the solutions with ^ — 7t/2 and ^ — 3tt/2 interchange 
the first and third columns, which does not agree with experiment. Thus, apart from a discrete ambiguity, the 
tribimaximal form is uniquely chosen when a smooth connection with the (23)-type perturbation is required. Of 
course, the smooth connection corresponds to choosing the correct initial states for the perturbation treatment. 

It is easy to see that perturbations which leave the other two S2 subgroups invariant, do not lead to mixing matrices 
of the desired tribimaximal form. The perturbation which commutes with S'^^^) is: 

t' u' O 

A' = ( u' t' O I . (19) 


Similarly, the perturbation which commutes with S"*^^"^^ has the form: 

/ t" O u" \ 

A" = O O O . (20) 
\ u" O t" / 

The stability condition for obtaining the tribimaximal mixing for the A' pertubation would require the matrix 
element {K^gJ^.JA'KTBM)l3 to vanish; instead it takes the value -^(f — 2u'). Similarly, the stability condition for 
the A" pertubation does not work since the matrix element {K^ g j^,^ A" Kt b m) 13 takes the generally non-zero value 
^{-t" + 2u"). 

While we have seen that the stability condition for (23) invariant perturbations enforces the experimentally plausible 
tribimaximal mixing, the underlying S3 symmctry should allow characteristic stable mixing matrices to emerge for 
either the (12) invariant or (13) invariant perturbations. What are their forms? In the case of a (12) perturbation, 
the stability condition associated with degenerate perturbation theory reads: 

{K^A'K)u = 0. (21) 

Here the characteristic mixing matrix emerges as K — KtbmX{^) for a suitable value of ^. The solution is easily 
seen to have the form: 



111 



KtbmX{-) = i ^ 71 v| i ■ (22) 

In the case of a (13) invariant perturbation, the stability condition associated with degenerate perturbation theory 
reads: 

-r 



(if^A"if)i3 = 0. (23) 



Here the characteristic stable mixing matrix turns out to be: 



-1 



,7r. 



V6 \/3 \/2 



KTBMXq) = I ^ ^ O I • (24) 



-1 

76 73 



The situation is summarized in TablelU Mathematically, any of the three perturbations will result in a stable mixing 
matrix. However, only the (23) perturbation gives the experimentally allowed tribimaximal form. For example, we see 
that the zero value of ifia, in good present agreement with experiment, only holds for the A [(23)-type] perturbation. 



III. ZEROTH ORDER SETUP 



In order to go further we adopt convenient conventions for the, in general, complex parameters a and (3 defined 
in Eq.([7|). The goal is to adjust a phase. V' in order that the zeroth order spectrum has three exactly degenerate 
neutrinos. As shown in Fig. [H we take the 2-vector 3/3 to be real positive. Then the 2-vector a lies in the third 
quadrant as: 

a = -^|a|e-*^/^ (25) 
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Perturbation 


Mixing matrix 


A 


Ktbm 


A' 




A" 





TABLE I: Characteristic, stable mixing matrices for each S2 invariant perturbation. 
where the physical phase ^ lies in the range: 

o < V < TT. (26) 

Finally |q;| is related to (3 by, 

I I 3/3 



' ' 2sm(V'/2)' 

In the limiting case ip ^ n, a takes the real value, 



(27) 




FIG. 1: Isosceles triangle with angle 1/) between the equal length 2-vectors a and a + 3/3. 



IV. ANALYSIS OF FAVORED PERTURBATION 



For simpHcity we will consider the parameters t and u in Eq. (|15p to be real rather than coinplex. The entire 
neutrino mass inatrix to first order is = al + /3d + A. Since we are working in a basis where the zeroth order 
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piece is diagonalized by the tribimaximal matrix, R, we must diagonalize the matrix: 

R^{al+ f3d + A)R = 
v 

3 



/ 



al + 



t + u ^{t + u) O 
^(t + u) 3/3+|(t + u) O 



(29) 



V 



o O t-u 

Diagonalizing the upper left 2x2 sub-matrix yields the three, in general, complex eigenvalues: 

3, ' 



:(/3 + r)(l 



a+-{p + T)il 
a + t ^ u, 



'l 



4/3T 



3(/3 + T)2 



4/3T 



3(/3 + T)2 



a + 3/3 + 2T, 



(30) 



where we introduced the abbreviation, T = (t + u)/3. The indicated approximations to the exact eigenvalues corre- 
spond to working to first order in the parameters t and u. Remember that according to our original setup, t and u 
are supposed to be small compared to \a\ and /3. Since Fig. [T]shows that generally \a\ > 3/3/2, it is sufhcient that \t\ 
and |m| be small compared to /3. 

In this approximation the corresponding eigenveetors are the columns of, 



1 ^{t + u) O 



1 O 
o 1 



(31) 



The entire diagonalization may be presented as, 



Here mi, m2 and ma are the three (positive) neutrino masses and 

K = RRiP 



771 1 











1712 











TO3 



(32) 



(33) 



is the full neutrino mixing matrix (in a basis where the charged leptons are diagonal). The neutrino masses, to order 
{t, u)/f3, are seen to be: 



3/3 



1- 77^(^ + w)sin^(-^) 



3/3 ^ 
3/3 



1 + -^{t + u)sin'^{^) 



1- —it-u)sin^{^) 



2 

4 
6 
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These mass parameters were made real, positive by the introduction of the phase matrix: 

/ 

P = 



where, 



(34) 




(35) 



TT 1 1, cot(ih/2) 



2 2 
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1 1 , cot{^l2) 



2 1 I 4(t+«) ^ 

"l' 9/3 

3/3 

To compare with experiment, we have important information from neutrino oscillation experiments (sj-p^. It is 
known that 1321 



B= |to2_j„2| ^ (2.5±0.5) X lO^^eF^. (37) 



A = - m{ = (8 ± 0.3) x lO^^eV^ 
2 2 1 

Also, constraints on cosmological structure formation yield [ssl j a rough bound, 

mi+m2+m3 < 0.7eV. (38) 

The two allowed spectrum types are: 

Typel : TO3 > m2 > mi, 

Type2 : TO2 > mi > TO3. (39) 

Now, from Eq. ([Ml) we see to leading order: 

ml - ml = 3l3{t + u), 

mj-mj = f3{-5t + u). (40) 
The quantities pt and Pu may thus be obtained for a type 1 spectrum as: 

(3t = A/m ~ B/6 w -4.13 X lO^-^e^^ 

/3u = 5^/18 + 5/6 w 4.39 X lO^^'eF^ (41) 

where the central experimental values were used. In the type 2 spectrum case, we should change B —* —B in the 
above to find, 

Pt = ^/18 + B/6 « 4.21 X lO-^eF^ 

/3u = 5/1/18 -B/6 w -3.94 X lO^^eV^ (42) 

Thus the, assumed real, violation parameters fit and fiu are now known for each spectrum type. Information 
about the quantity /? may in principle be obtained from the perturbed lepton mixing matrix given in Eq. (j33p : 



K . 



/ -2 _ \/2(f+n) 2(t+M) Q \ 

\/6 9/3V3 V3 9/3\/3 

\/2(t+M) 1 , (t+n) 1 

76 9/3^3 73 9/3^3 72 

\/2(f+M) 1 I (t+ii) -1 



P. (43) 



\ 76 9/373 73 9/373 72 / 

With a usual parameterization (3^ the matrix with zero (13) element takes the form, 



Cl2 Sl2 O 

K = I -S12C23 C12C23 S23 I P, (44) 



S12S23 -C12S23 C23 
where C12 is short for cosdn for example. This amounts to the predictions, 

2 ^/2{^it + /3u) 



C12 



V6 9v^/32 ' 



si3 = 0. (45) 
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Notice that, when the perturbation is absent, this agrees with the tribimaximal form used here if both 612 and 6*23 lie 
in the second quadrant. The results of a recent study fs^ of neutrino oscillation experiments are: 

{s,,f = 0.304tS]°?^. (46) 

One immediately notices that the prediction, (523)^=1/2 is unchanged from its tribimaxinial value by the perturbation 
and agrees with the new analysis. On the other hand the tribimaximal prediction, (si2)^=l/3 is sHghtly changed 
from its tribimaxinial value and actuaUy hes shghtly above the upper experimental error bar. This is probably not a 
serious disagreement but it might be instructive to try to fix it using the predicted perturbation in the present model: 

(^12) -3-^^^- (47) 

For either the type 1 or type 2 assumed spectrum, the perturbation is seen to be in the correct direetion to lower the 
value of 512)^, as desired. However, because of the large cancellation between (it and /3u, this effect is extremely small 
for a reasonable value of /3^; even with (3 as small as 0.05 eV, (512)^ is only lowered to 0.332. 

It is also interesting to discuss the absolute masses of the neutrinos rather than just the differences of their squares. 
Since the differences are known, let us focus on one of them, say 7713: 

"^3 ~ -i^csc{-) sin[-). (48) 

Notice that the first term on the right hand side is, using Eg.ipT]). simply the zeroth order degenerate mass, |a| while 
the second term represents the correction. Also note that (see Fig. [T]) the point V = O is not allowed. Considering 
■0 as a parameter (related to the strength of neutrinoless double beta decay), this equation represents a quadratic 
formula giving (3 in terms of the absolute mass for any assumed ip. In Fig. [21 adopting the criterion that \t\/ (3 and 
\u\/ [3 be less than 1/5 for perturbative behavior, we display the perturbative region in the ms — ?/; plane. In contrast 
to the case of 777.3, "^i and 7772 are seen to have small corrections since they of course depend on /3(i + u) rather than 
I3{t~u). 



V. NEUTRINOLESS DOUBLE BETA DECAY 



The characteristic physical novelty of of a theory with Majorana type neutrinos is the prediction of a small, but 
non-zero, rate for the neutrinoless double beta decay of a nucleus: (A, Z) — > (A,Z+2) + 2e^. The appropriate leptonic 
factor describing the amplitude for this process is, 

|777ee| = \mi{Kllf + 7772(^^12)' + 7773 (if 13)' | • (49) 

Substituting in the neutrino masses to order (t, u)/ 13 from Eq . ([Ml as well as Ea. ([55)) yields: 

,™ 1^ 3/? ,2 , 4(t + 77) 2>.^.l 4(^ + 77) 2.V'^i 2^(r-.^| /.r,. 

cos (-) + [- Ti^^cos {-)]e ^ )\. (50) 



' ' 2sm(|)'3 27/3 ^2^ ^3 27/3 

This expression can be straightforwardly evaluated. The needed intermediate quantity cos[2[t — cr)]may be easily 
obtained from Eqs. ([36)) with the construction of a suitablc right triangle. We find, correct to first order in {t,u)/f3: 

I777..I « f C.c(f )y5TW[l + ^^^^]- (51) 

We recognize the first term as the ^3 limit result, Eg. p^ . Note that due to the large cancellation in the quantity 
{t + u), the first order correction to |777ee| is strongly suppressed. 
The experimental bound on |777ee| is given [s^ as. 



|777ee| < (0.35- 1.30)ey, (52) 

which is small enough so that there is hope the possibility of a Majorana neutrino might be settled in the near future. 
Since the correction to |r77ee| has been seen to be small in this model we can take over the zeroth order inequality 
in Ea. fTlj) . This means that the existence of the Majorana phase, ip can alter the amplitude for neutrinoless double 
beta decay by a factor of three for given (approximately degenerate) neutrino masses. 
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FIG. 2: Sketch of perturbative region in the m^-il) plane. It is about the same for both type 1 and type 2 neutrino spectra. 
Note that is measured in radians and mz is measured in eV. 



VI. SUMMARY AND DISCUSSION 



In some ways the problem of "fiavor" in the Standard Model is reminiscent of that in Strong Interaction physics 
before the quark model. At that time it was realized that, as a precm'sor to detailed dynamics, group theory might 
give important clues. 

Then the strong interactions were postulated to be SU(3) fiavor invariant with a weaker piece having just the the 
SU(2) isospin (times hypercharge) invariance. In addition it was known that there was a still weaker isospin breaking 
(possibly QED) which by itself preserved a different SU(2) invariance (so-called U-spin). 

Here, an analogy for neutrinos of the first two steps was studied in a perturbative framework. In T| and [2] the 
possibility that the neutrinos are not strictly degnerate at the unperturbed level was contemplated. However in this 
paper we have examined a strictly degenerate first stage (setting to zero the parameters called respectively e and b in 
those two papers). 

At the 5*3 invariant level the neutrino mixing matrix is actually arbitrary up to a rotation in a 2-dimensional 
subspace. This problem can be settled (since degenerate perturbation theory is involved) by specifying the transfor- 
mation property of the perturbation to be added. Although there is widespread agreement that the first perturbation 
should preserve the ^2 subgroup which interchanges the sccond and third neutrinos, we presented in section II, for 
completeness and interest, the mixing matrices for the other two possibilities also. 

In sections III and IV we carried out the perturbation analysis for any choice of a Majorana-type phase, i/" which 
plays an important role in this model. If il) is considered fixed there are three parameters in the model (In |2(] ij} was 
considered fixed at the value tt.) These three parameters can be taken as (3t, f3u and P defined above. The quantities 
f3t and pu were found in terms of the neutrino squared mass differences for each choice of neutrino spectrum type, 
ie normal or inverted hierarchy. The value of /3 depends on the presently unknown absolute value of any neutrino 
mass. The magnitudes of pt and f3u are similar (though not exactly equal) but differ in sign. Thus the perturbation 
corrections which involve (/3t + /3u) are very small. Clearly (see the first of Egs. ipO)) ) this is due to the small solar 
neutrino mass difference. This situation occurs for the correction to the mixing parameter sin^9i2 in addition to mi 
and 7712, the masses of the first two neutrinos. The perturbation dependence on -Pu) is not suppressed however. 
This occurs for the mass, of the third neutrino. This result was used to make a sketch of the region in the ^p-m^ 
plane for which the perturbation approach given seems numcrically reasonable. 
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In section V, the perturbation was calculated for the neutrinoless double beta decay factor, niee, going beyond 
the formula obtained at zeroth order in [l[. The correction turns out to be proportional to {(3t + (3u) and hence is 
suppressed. However the effect of the Majorana phase in this model can be rather important. 

The explicit role of the Higgs sector, which is believed to be at the heart of the matter, was not discussed in the 
present paper. However, this as well as some further technical details were discussed in [Ij]. Further treatment of 
this aspect is interesting for future work as is a detailed investigation of the weakest perturbation, the analog of the 
U-spin preserving perturbation in the strong interaction. This could be used to further study other consequences 
of possibly non-zero ^13. 
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